Moyal phase-space analysis of nonlinear optical Kerr media 



T. A. Osborn 1 and Karl-Peter Marzlin 2,3 

'Department of Physics and Astronomy, University of Manitoba, Winnipeg, Manitoba R3T 2N2, Canada 
2 Department of Physics, St. Francis Xavier University, Antigonish, Nova Scotia, B2G 2W5, Canada 
3 Department of Physics and Astronomy, University of Calgary, Calgary, Alberta T2N 1N4, Canada 

Nonlinear optical media of Kerr type are described by a particular version of an anharmonic 
quantum harmonic oscillator. The dynamics of this system can be described using the Moyal 
equations of motion, which correspond to a quantum phase space representation of the Heisenberg 
equations of motion. For the Kerr system we derive exact solutions of the Moyal equations for a 
complete set of observables formed from the photon creation and annihilation operators. These 
Moyal solutions incorporate the asymptotics of the classical limit in a simple explicit form. An 
unusual feature of these solutions is that they exhibit periodic singularities in the time variable. 
These singularities are removed by the phase space averaging required to construct the expectation 
value for an arbitrary initial state. Nevertheless, for strongly number-squeezed initial states the 
effects of the singularity remain observable. 



I. INTRODUCTION 



Phase space methods provide one of the most important tools to investigate the relation between classical and 
quantum mechanics. In classical dynamics, the state of a system can be described by a probability distribution that 
is a function of x = (q,p), i.e., of position q and momentum p of a particle. In quantum physics such a distribution 
cannot exist because Heisenberg's uncertainty relation prohibits simultaneous knowledge of q and p, but a number of 
quasi-distributions have been proposed for a phase space analysis of quantum systems p] . Among the most popular 
is the Wigner function which is related to the density matrix p(q,q') = {q\p\q') of a particle by 

1 f°° 

W(x) = — / p(q-q',q+ q') e 2 ™ /« dq' . (1) 

Usually the parameter £ is replaced by h, but for reasons explained below we use £ instead. The Wigner function is 
not a probability distribution because it can have negative values, which are often interpreted as an indication for 
genuine quantum effects. 

The Wigner function and its time evolution for a given Hamiltonian have been applied to analyze a huge variety 
of phenomena; in particular during the last two decades it has been used to analyze the quantum state of light [2] . 
However, to investigate the classical-quantum correspondence the closely related concept of the Weyl symbol [A] w of 
an operator A may be more suitable. The Weyl symbol has essentially the same definition as the Wigner function 

t, with p replaced by a general operator A. For the special case of the density matrix it leads to the relation 
w = 2tt^W{x). The difference between the Weyl symbol and the Wigner function is that the latter refers to the 
quantum state while the former refers to observables. For an observable, the time evolution of its Weyl symbol is 
governed by the Moyal equation 



d 



-A w (t) = {A w (t),H w } M . (2) 



Above H w ,A w (t) are the Weyl symbols of the Hamiltonian and A(t), respectively. The Moyal bracket {-,-}m, cf. 



(A4), is the quantum extension of the classical Poisson bracket. The Moyal equation corresponds to a phase space 
formulation of the operator valued Heisenberg equations of motion. The similarity between the Heisenberg equation 
for observables and the classical equations of motion makes the Weyl symbol representation a powerful tool to shed 
light on the relation between classical and quantum dynamics. The Moyal equation of motion enables one to study 
the dynamics of observables without reference to the quantum state of the system. 

Despite these advantages the Moyal equation has not been used extensively to analyze quantum systems because 
it is considerably more difficult to solve than the Schrodinger equation. For this reason not many exact solutions 
are known, and a comparison with experimental data is often difficult. The purpose of this paper is to improve this 
situation by providing an exact solution of the Moyal equation for an experimentally relevant system: the Kerr model 
of nonlinear optics for a single-mode of the quantized radiation field. We obtain analytical solutions of the Moyal 
equations that are valid for all time. These exact solutions are then used to characterize the transition from classical 
to quantum dynamics and to provide a phase-space based physical interpretation of the Kerr effect. 
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This paper is organized as follows. In Sec. [TT] we review the basic properties of the Kerr effect. A closed, exact 
expression ( 16 1 for the phase space representation of all Hciscnberg operators in the quantized Kerr model is derived 
in Sec. |III| and its classical limit is studied in Sec. |IV| In Sec. [V] we show that the periodic divergences in the 
Moyal representation of the annihilation operator generate characteristic, finite peaks in the expectation values of 
the canonical variables for squeezed coherent states of light. In Sec. VI it is shown that the expectation values of all 
observables remain finite for general initial states. 



II. THE KERR MODEL 



The Kerr model of optical nonlinearities is one of the most studied systems in quantum optics. In a Kerr medium, 
the refractive index of a classical beam of light depends on the light intensity as n = + n^I, where n denotes 
the total refractive index, the linear refractive index, / the light intensity, and the optical Kerr coefficient. 
The Kerr effect is invaluable for spectral broadening and self- focusing of laser pulses [3]. It usually appears in 
special crystals and its magnitude is typically so small that large light intensities are needed. However, recent 
research on electromagnetically induced transparency [H |5J |5] has made very large Kerr coefficients with values 
of up to w 0.1cm 2 /W possible [3 [H [U [TU] and may even lead to nonlinear effects at the single-photon level 

PHH2 031 HUE]. 

A quantum description of the Kerr effect can be accomplished by replacing the intensity of light in the refractive 
index by the corresponding operator. This is equivalent to introducing a quartic interaction term in the radiation 
Hamiltonian £Q . A particularly simple description can be achieved if the photon dynamics is confined by an optical 
cavity with high finesse mirrors. Such cavities may support only a single light mode in a given spectral range, so that 
the dynamics can be described by operators a, that annihilate or create a photon in the cavity mode, respectively. 
The photon number operator is given by N = a* a, and the Kerr Hamiltonian is given by the Wick ordered operator 

H = w 2 (a t ) 2 a 2 +cjia f a . (3) 

Physically, lo\ is related to the linear index of refraction by u>% — kv gI n^ and L02 to the nonlinear refractive index by 
lo 2 = kv gr n^ I Q , where w gr is the group velocity of light in the medium, k its wave number, and Iq = 2hkc 2 /V the 
intensity of a single photon in a cavity of volume V. The operators H,a act on the Hilbert space H — L 2 (U.,C) and 
satisfy harmonic oscillator commutation relations 

[a, a f ] = £1 [N, a] = [N, a f ] = £a f . (4) 

With this notation we have introduced the real dimensionless parameter £ that allows us to interpolate between 
classical and quantum evolution. The fundamental distinction between quantum and classical mechanics resides 
in commutivity. The product of observables in classical mechanics is abelian whereas the product operation in 
quantum mechanics is noncommutative. As £ — * 0, noncommuting behavior in the Kerr model is suppressed, and for 
£ = 1 standard single mode photon physics is recovered. Specifically, the Moyal equation of motion automatically 
incorporates Bohr's correspondence principle: in the limit £ — > 0, the Moyal bracket becomes the Poisson bracket and 
Eq. Q then turns into the Poisson equation of motion. 

The relation between £ and the conventional "quantization parameter" h can be seen by relating the creation and 
annihilation operator to two Hermitian operators via a = (q + ip)/y2 and — (q — ip)/y2. For dimensionless 
position and momentum operators x = (q,p) = $2), one then has 

[xj,x k ] = i£J jk J ^ -°i ) ' ( 5 ) 

where J is the Poisson matrix. In a coordinate parametrization where q is proportional to length and p to momentum, 
Xk] = ihJjk, so the limit £ — > is equivalent to letting Planck's constant h go to zero for a conventional quantum 
harmonic oscillator. We use £ instead of h because in quantum optics the operators a, a' have a different physical 
interpretation than for a single Schrodinger particle in a harmonic potential. As a consequence, their degree of 
commutativity is not controlled by h, but rather by the mathematically introduced deformation parameter £. 

Throughout our derivations free use is made of the Weyl symbol calculus that represents Hilbert space operators 
by functions in phase space. An overview of this quantum phase space representation and its non-commutative * 
product is presented in the Appendix A. 
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III. THE MOYAL-KERR PROBLEM AND ITS SOLUTION 



In this section we describe Heisenberg picture evolution in Weyl symbol form, identify the symmetries of the Moyal 
equation of motion and use these symmetries to construct an exact solution. 
The Weyl symbol representation of the Haniiltonian ^ is 



[H] w = H(S,x)=W2 



(6) 



where x 2 = q 2 + p 2 . The u)2 term is the phase space form of the non-linear interaction. The uj\ portion is the symbol 
of the number operator [iV] w (a;) = N(x) — \{x 2 — £) and represents the evolution of non-interacting photons. 



First consider the Moyal equation for a general observable. Denote Schrodinger evolution by Ut 
Let Oq be an observable with dynamical value Q(t) = U} Oq Ut and Heisenberg equation 



exp(-itH/£). 



The Weyl symbol image of Eq. Q is Moyal's equation Let Q(t\x) 
observable, then 



(7) 



[0(t)] w (a;) be the symbol of the evolving 



e(t\x) = {e(t),H} M (x)^ic 1 (H*e>(t)~e(t)*H)(x) 

= iC 1 (H(L) - H(K)) e(t\x) . 



(8) 



In the last identity, we employ the expression of the Moyal bracket in terms of the left and right operators L and TZ 
which are defined in Eq. ( A 7 ) . It converts the Moyal bracket into a differential operator acting on the target function 
G(t\x). Evaluating, H(L) - H(TZ) for Haniiltonian ^ one obtains the following third order differential equation 



6(t|a;) 



uj 2 x 



2£ 



f 2 

^d 2 
4 x 



[x ■ Jd x )G(t\x) 



(9) 



To fully characterize a quantum system, Eq. ([9| has to be solved for a complete set of operators. For the Kerr-Moyal 
problem, such a set is given by {(a') s d m \0 < s,m € N}. We denote by 



e 8m (t|a:) = [{a{tf) s a{t) m ] w {x) 



(10) 



the Weyl symbol of the corresponding operators in Heisenberg picture, where a(t) = U t aUt- Our task is to solve 
Eq. ^ for the set of symbols Q sm (t\x) with initial conditions Q sm (0\x) = [(a^(0)) s d m (0)] w (a;). 

Moyal equation |9]) has the form of a Schrodinger equation over the cc-variable manifold. Specifically, the function 
®sm(t\x) may be considered an unnormalized 'state' over the manifold T*R = M 2 . This is a general feature of the 
Moyal equation and has been used to construct WKB type asymptotic approximations [TB] for Q(t\x). The system ^ 
admits a standard [17 small £ expansion because the highest order differential operator, the Laplacian d 2 , is scaled 

Quantities like L and TZ act on the Weyl symbols <d sm (t\x), while operators like a act on the usual Hilbert space. 
To distinguish between these two cases we use a hat to denote the latter and script capital letters denote operators 
acting on the Hilbert space Ti.2 = L 2 (M 2 , d 2 x) . 

It is now useful to determine the symmetries present in the equation of motion ([9| . In this context it is advantageous 
to introduce complex coordinates z = (q + ip) and d z = \{d q — id p ). In this notation [a] w (x) = a(x) — zj\/2 and 



[d^] w (x) = a(xY 



z*/-\/2. Because a is a linear function: a-ka = a 2 , etc., giving 
[a m ] w (x) = 2- m ' 2 z m , [(a^YUx) = 2~ s ' 2 z*' 



Employing Eq. (A8l it is then straightforward to show that the initial condition is 

1 



e sm (0\x) = d(L) s a m (x) = 



min(s,m) 



V2 



X W{m,s,l)(-l) l a s -\x)a m -\x) 



(11) 



1=0 
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with W(m, s, I) = s\m\ [l\(s — l)\(m — /)!] 1 . For the set of operators under consideration, the Kerr-Moyal equation 
|9]) then becomes 

e sm (t\x) = - [u 2 K + wi] {x ■ Jd x )Q sm {t\x) (12) 

with 



[x -Jd x ) = i(zd z - z*<9 z .) . 

The action of the phase space operator (x • Jd x ) is similar to that of the angular momentum operator L z on Hilbert 
space. This can be easily seen in polar coordinates, z — re 1 ^, where it takes the form [x ■ Jd x ) — . Furthermore, 
the initial condition Q sm (0\x) is an eigenstate of (x ■ Jd x ) with eigenvalue A sm = i(m — s). Because [fC, (x ■ Jd x )] = 
we can infer that O sm (t|x) will remain an eigenstate of (x • Jd x ) with the same eigenvalue. Via this eigenfunction 
mechanism the third order partial differential equation (12) is reduced to second order, 



®sm{t\ x ) = -i(m - s) [uj 2 fc + ^1] <3>sm(t\x) . 



(13) 

We remark that if m = s the right side of Eq. (13) is zero. This means that Q mm (t\x) — (x 2 /2) m is a constant of 
motion. This function is also a classical constant of motion because {H, x^ 2 ™-*} = 0. 
Equation (13 1 has the formal solution 



e sm (t\x) = e- i{m - s)uit e- l( - m - s)tul2K <d sm {Q\x) 



We can now take advantage of the special form ( 11 ) of the initial conditions. It is well known cf. 



(14) 

, pg. 40) that 



A = e r Ae K for [K, A] = rA and r e C. Because 



we can express the formal solution as 

e sm (t\x) = 2- s / 2 exp(~i(m- S )t(L0 1 + Lj 2 S; S ))(z* -Zd z ) s e-^ m - s ^ K f-=J 

= 2~ s/2 exp ( - i(m - s)tLu 2 £s) (z* - £d z ) s e 0m (m- 1 (m - s)t\x) . 
It is therefore sufficient to find a closed form for <do m (t\x). To do so we make the ansatz 

Q 0m (t\x) = e 9{t)x \- lm{ ^- 2 ^ ]t f{t)a m {x) , (15) 

with initial conditions g(0) = and /(0) = 1. Inserting this ansatz into Eq. ( |13[ ) and sorting the resulting equation 
in powers of \z\ 2 yields a coupled set of differential equations for g(t) and f(t), 

g = imu} 2 (-l + £ 2 g 2 ) f = im(m + l)£ 2 uj 2 gf , 



which have the solutions 



<?(*) 



tan(m^o;2^) f(t) — (sec(m^u; 2 t)) 



m+l 



This leads to one of the main results of this work: the exact solution of the Moyal equation for the complete set of 
operators (a^) s a m is given by 



Q sm {t\x) 



-i(m-s)iu 1 t e i(2-s)t/ gec £\ro+l 



V2 



exp 



— zz tani 



(16) 



with t — (m — s)^uj2t- Evaluating the s-fold derivative and converting to the phase space variables x yields 



i(s,m) 



e sm (t\x) = e-^ m - s ^ lt (seci) s+m+1 exp (2U- tant) J] W(m,s,l) 



cost) 1 a s - l (x) a rn - l (x) . (17) 



1=0 



This final form displays the adjoint symmetry: [(d(<:)^) s a(i) m ]^ — (a(i)^) m a(i) s , or equivalently Q sm (t\x)* = Q ms (t\x). 



A striking feature of the solutions (17 1 is that they have a singular amplitude for times whenever cost = 0. 



Henceforth we will refer to this behavior as the Moyal singularity. Its mathematical origin is that g(t) obeys a 
non-linear Ricatti equation. 
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IV. CLASSICAL AND QUANTUM TRAJECTORIES 

The manner in which quantum phase space solutions embed the classical dynamics occurs in two different ways. In 
the first way, one characterizes how the solutions <d sm (t\x) transform into the Poisson equation solutions as £ — > 0. 
The second semiclassical association relates quantum expectation values to corresponding classical flows. The first 
way, the phase space correspondence, is treated in this section. 

Quantum trajectories on phase space are defined as the symbol image of the Heisenberg coordinate operator 
evolution: x(t) = llJxUt, in detail 

Z(t, eias) = [x(t)] w (x) = ([g(i)] w , W)U) 0*0 ■ (18) 
The Moyal solutions above give formulas for Z(t,£\x) via the relationships 

[q(t)Ux) = ~(e 01 (t\x) + G w {t\x)) = V2Ree i{t\x) 

\p(t)] w (x) = ^=(e 01 (t\x) - e 10 {t\x)) = V2lmQ 01 (t\x) . 

The dynamics for a the classical version of the Kerr problem is simple. The £ = part of H(£, x) defines the 
classical Hamiltonian 

#(£, x) = H cl (x) + £hi(x) + ^h 2 (x) 
H c \{x) = \oj 2 x a + \u\x 2 hi(x) = ~lu 2 x 2 h 2 (x) = \w 2 ■ 

The classical trajectory Z c \(t\x) = (^ c i l^) ^ Pci(^ |^)) is then the solution of Hamilton's equation 

Z cl (t\x) = Jd x H cl (Z cl (t\x)) 
with initial condition Z c i(0|ar) = x. The solution is Z c \(t\x) — [exp t{uj 2 x 2 + u>i)J] x. In matrix form this is 

Zoi(t\x) = ( q *ffi\ ) = ( «?(^ + "0 * M"2X 2 + "l)t \fq\ (1Q) 
\Pcl(t\x) J y—Sm[UJ2X +Wi)t C0S(W22T + U>i) t J \ P J 

This is oscillatory motion with a variable frequency u 2 x 2 + u)\ that depends on the initial value, x. The frequency 
increases as the constant of motion x 2 increases. 

To compare the classical and the quantum trajectory it is useful to introduce the complex quantity 

a cl (t\x) = ^=(q cl (t\x)+i Pcl (t\x)) = e -^ 2 +- 1 )*-L(g + ^), (20) 

which is the classical quantity corresponding to the annihilation operator. It agrees with the predictions of the Kerr 
model for the complex electric field amplitude of classical light: the leading phase factor represents the phase shift 
that light would experience when it travels through a nonlinear medium of length L = tv gr . In an optical system, x 2 
represents the mean number of photons in the cavity, which can also be expressed as the intensity of the light field in 
units of the intensity of a single photon in the cavity. 

Employing Eq. (19 1 we can express the related quantum trajectory as 

e„i(t|a;) = sec 2 (^ 2 i) e^^a cl (t\x) , (21) 

with the quantum phase factor 

x, t) = 2& 2 t + x 2 {u 2 t - C 1 tan(f w 2 *)) . (22) 



This phase vanishes at £ — 0. If one implements a power series expansion of (21 1 about £ = 0, the result defines the 
semiclassical expansion of the Moyal solution. To first order in £ one has 

e s c 1 (t\x) = a cl (t\x)[l + 2iLU 2 t£ + O(e)] ■ (23) 

Expressions ( 2l|23 1 show how the periodic quantum and classical flows are interdependent. If \£,co 2 t\ -C 7r/2 and 
|3>(£, x,t)\ <C 7r/2 the quantum and classical trajectories nearly coincide. This limit very well describes all experiments 
with conventional nonlinear optical crystals for which the nonlinear refractive index is very small. 



FIG. 1: Magnitude of the ratio between the quantum amplitude Oqi and the classical amplitude a c i. 




FIG. 2: Quantum phase factor <&(t;,x,t) for £ = 1. 



At the other extreme, when £,u>2t is close to an odd multiple of 7r/2, the amplitude factor sec 2 (^u>2t) is diverging and 
the quantum phase rotation x, t) is undergoing near infinite oscillation. This regime should soon be experimentally 
accessible by using EIT-based nonlinear media El HE HQ] , but below we will explore if such large oscillations are 
actually observable. 

Fig. [I] shows the quantum effect of the nonlinearity on the complex field amplitude |oci(i|x)|, which is independent 
of x 2 . It can clearly be seen that the periodic divergences disappear in the classical limit £ — > 0, and that even for a 
fully quantized theory (£ = 1) they only appear for very large nonlinear refractive indices. Fig.[2]displays the quantum 
phase factor <f> for a fully quantized theory (£ = 1). Like the amplitude it displays a periodic divergence in time. The 
width of the divergences in phase space is proportional to x 2 , indicating that they are an intensity-dependent effect. 
The divergences disappear in the classical limit £ — > 0. 

Generally, not everything is rapidly oscillating. Recall that x 2 is both a classical and quantum constant of motion. 
At the classical level this constant is recovered from the flow via Z c i(t, £,\x) 2 — x 2 . At th e qu antum level one has 



Z(t,£) * Z(t,£)(x) — x 2 . This latter identity can be derived from Berezin's representation (A5l of the star product 
which leads to a four dimensional Fresnel integral with value x 2 . 

In the literature |19l 120] the quantum trajectory Z(t,£\x) is generally approximated by a small £ asymptotic 
approximation 

Z{t,f\x) = Z cl (t\x) + tzV>(t\x) + ^ 2) (t\x) + ■■■ (24) 

In the Kerr-Moyal problem, and for most problems, this expansion has terms to all order in £. The exception occurs 
when H is quadratic. In this special case just the leading term Z c \(t\x) is nonzero. For this reason little can be learned 
about the general nature of the classical-quantum transition by investigating quadratic Hamiltonian problems. This 
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feature is seen in Moyal solutions above. The effect of the ^-quadratic part of Hamiltonian (|6| on Q sm (t\x) is confined 
to the phase factor exp(— i(m — s)LJ±t). This factor has no £ or x dependence. 

If an exact formula for Z(t,^\x) is known then these higher order expansion coefficients are given by 



z^{t\x) = ^Z{U\x) 



For example, the leading semiclassical correction to the classical Kerr problem flow is 

d 



(t\x) 



-2u! 2 t Z ci (t\x) 



(25) 



The method for computing z^ n '(t\x) when the full quantum trajectory is not available is to expand the Moyal 
equation identity in powers of £. This approach works if both the observable and the Wcyl system Hamiltonian are 
semiclassical admissible, namely both admit a power series expansion about £ = 0. This is the situation in the Kerr 
problem. The 0(£ 1 ) portion of the Moyal equation for Z(t,£\x) is an inhomogeneous Jacobi field equation for the 
unknown z^(t\x), i. e. 



Jt -JH'a(Z cX (t\x)) 



&\t\x) = Jd x h 1 {Z*{t\x)) 



(26) 



with initial condition z^^Ola;) = 0. The quantity H c is the Hessian matrix of H c . Simil ar e quations define the higher 
order corrections z^ n \t\x). A Jacobi field is a solution the homogenous version of Eq. (26 1 and provides a linearized 
prediction for small deviations about the classical flow Z c \(t\x). 

One can readily check that z^(t\x) — —2uj2tZ c \(t\x) is a solution of Eq. (26). This demonstrates the compatibility 

(IM 



of Z(t, Z\x) in Eq. (18 1 with the standard asymptotic semiclassical expansion generated by Eq. p4| . Formula (25 1 
also illustrates the small time limitation of this expansion. The z^ correction has unbounded growth in t; in order 
for the correction to be small one requires |2w2*| *C 1. The next correction term z^ grows like x 2 (u)2t) 3 . This shows 



that the expansion (24 1 is non-uniform in the (t,x) G 



domain. 



V. DYNAMICAL EXPECTATION VALUES AND THEIR CLASSICAL LIMIT 

In this section we compute the squeezed state expectation value of the Moyal solution corresponding to q(t) and p(t) 
and characterize their semiclassical limits. The squeezed states are of particular interest because they will allow us to 
study the effects of the singularity. Furthermore, squeezed states are of high practical value because they correspond 
to non-classical states of light which can be used for quantum information |21j . In addition, the reduced noise of 
specific observables makes them of interest in high-precision experiments such as gravitational wave interferometers 

m\ 

Squeezed states are unitary modifications of coherent states. We recall the defining equations for coherent and 
squeezed states. The coherent states \a), are translated vacuum states (see, e.g., Ref. [T]). The translation operator 
D(a), aeC, arga 6 [0, 2ir) shifts a by 

(a)aD(a) = a + al , D(a) = exp [£ -1 (aa^ — a*a)] . 

Defining \a) — D(a)\0), it follows that a\a) — a\a). The states \a) have unit normalization with inner product 

(a|/3) = exp[-r 1 (|H 2 + ||/3| 2 -a*/3)] . 

Given the coherent states one obtains the squeezed states by the action of a unitary Bogoliubov operator, V(t) = 
exp[r(a' t ) 2 — r*a 2 ], r = |r| exp(z</>), <j> £ [0, 2ir). The squeezed states \ra) = V(r)\a) are the eigenfunctions of 



V(t) aV( T y 



ra) = a\Ta) 



Because V(t) is unitary, the coherent state eigenvalue a and normalization (tq;|t/3) = {a\(3) are unchanged. 

The coherent and squeezed states are interpreted as being near classical because they have special properties with 
respect to the uncertainty relations. The squeezed state mean values are readily found by exploiting the metaplectic 
nature of the Bogoliubov transform, specifically, 

/ s cos 2 (6/2) + s- 1 sin 2 U/2) -I (s" 1 - s) sm(6) \ 

V(r)xV(r)^S(r)x S(r) = V^-s)^) ^ -s 2 (0/2) + s sin 2 (6/2) ' ^ 
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Im(a) 



arg a = n/2 
Aij) =0 



Re(a) 



arg a = 
A(f> = ji 

m - Jt 



FIG. 3: Representation of squeezed states in the plane of complex amplitudes. Each ellipse is centered at the mean value (a) 
with the principal axes corresponding to the uncertainties of the canonical variables. The length of the axes is determined by 
the squeezing parameter s and the direction of the ellipse by the phase factor (j>. The parameter A</> appears in Eq. (351. 



The parameter s = e xp(— 2£|t|) < 1 describes the amount by which the uncertainty of a canonical variable can be 
reduced (see Eq. (28 1). The set {S(t)\t S C} is a family of positive symplectic matrices with inverse5'(T)" 1 = S(—t) 
and group multiplication law S(t) 2 = S(2t). 



The q, p variances turn out to be |23j 



t ( 1 



■ cos 2 (0/2) + s 2 sin 2 (0/2) 



1 



(Ap) 2 Ta = £ S 2 cos 2 (0/2) + -sin 2 (0/2) 



The uncertainty statement appropriate for this context is the Schrodinger-Robertson inequality: 

e 



(Fy 



F = {q-(q),p-(p)} s 



(28) 



(29) 



with the anti-commutator {X, Y} sym = XY + YX. Employing (27 1 to evaluate F gives (F) Ta = (£/4) (s~ 2 — s 2 ) sin <fi. 

Combining these statements shows that the rev squeezed states are minimum uncertainty states with respect to the 
Schrodinger-Robertson lower bound. In fact [24J a state that fulfills the equality in (29 1 must be a squeezed state. 
We remark that the phase 4> of the squeezing parameter r determines which of the canonical variables is squeezed: 
for <f) — (4> = 7r) the variance of p (q) is reduced by a factor of s, respectively. In other words, the angle variable </>/2 
rotates the semi-axis of the uncertainty ellipse with respect to the q,p axis. This behaviour is usually visualized by 
representing the squeezed state as an ellipse in the complex a-plane that indicates the uncertainties of the canonical 
variables, see Fig. [3] 



The quantum phase space representation of the expectation value is the phase space integral ( A9 1 

1 



{a{t)) Ta = (Ta\a(t)\ra) = 



2tt£ 



Q m (t\x)[\Ta){Ta\] w (x)d 2 x 



where 0oi(t|x) is the symbol of a(t) = l/j aU t . 

Next we compute the Weyl symbol of [|ra)(ra|] by relating it to the simpler quantity [|a)(a|] i 
coherent state wave function 



(30) 



From the 



1/4 



exp 



2 



a q — a Re a 



(31) 



one obtains the associated Wigner distribution. Let x = (q,p) = (y/2~Rea,y/21ma) be the a coherent state mean 
values, then one has 



[\a)(a\]jq,p) = 2exp£ _1 {-(q 2 + p 2 ) + 2qq + 2pp-(q 2 +p 2 )} 



(32) 



The Weyl symbol (32 1 is real because |a)(a| is hermitian. The squeezed generalization of this follows from the Weyl 
symbol covariance property (AlOl 



[\ra}(ra\]Jx) = [V (r)\a) (a\V (r)^] Jx) = [\a)(a\] JS(t) x 
1 



2 exp 



5(2r) x + 2x ■ S(t) x-x- 



(33) 
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The density matrix |ra)(ra| is projection operator that characterizes a pure ensemble of photons with mean number 



[a)a) Ta = sinh^f |t|) + \a cosh(2£|r|) + aV sinh(2£|r|) | * . 

The photon number is a constant of motion since a* a commutes with the Kerr Hamiltonian. 

The integral (30 1 is conveniently computed by diagonalizing S(t) and S(2t). The phase space rotation 



R{4>) 



achieves this via 



S{t) = R(<p)A(s)R(- 



A( S ) 



s 
s" 1 



Note the S(t) eigenvalues Ai = s, A 2 = s^ 1 are independent of cj) and likewise the matrix R(<fi) is independent of s. 
The eigenvalues are positive because S'(r) > and are mutual inverses since det5'(r) = 1. 
Introduce the variables y = R(—<j>) x and 



V = R(-<i 



_ ( gcos(0/2)+psrru>/2) \ _ ( q<p 
-q sm(<j>/2) +p cos(0/2) 



Changing the integration variable from x to y and employing the identity ®ai(t\R(4>)y) — exp(i0/2) Ooi(t|y) allows 
one to write the integral ( 30 1 as 



{a(t)) Ta 



1 



9 i(t|y)exp9 - y K{s 2 )y + 2y A{s)y-y 2 



(34) 



In displaying the final result it is useful to use the abbreviations 



T = tan(£w 2 t) , G(T, s) 



(1 + iT)- 



(l + is- 2 T)(l + is 2 T) 



Integral ( 34 1 is a generalized Gaussian integral and evaluates to 



(a(t)) Ta =aG 3 / 2 



- cos f£w 2 i 
s V 



-— ) +t8B m[^t-—) 



exp {-i^t +i^t- A(j)/2)) 



xexp ^~ 2l i(ITW 



iT + s~ 2 cos 2 (A</)/2) + s 2 sin 2 (A0/2) 



(35) 



Here A0 = — 2 arg(a) is the difference between the squeezing angle and twice the phase of the coherent state 
amplitude a [35]. The branch cut for vG lies along the positive real axis. 

The result above constructs the quantum mean (a(t)) Ta directly from the Moyal solution Q 01 (t\x). It describes in 
detail the dependence of the expectation value on the semiclassical scaling parameter £ as well has the squeezing and 
coherent state variables, r and a. The (q(t)) Ta and (p(t)) Ta predictions are obtained from the real and imaginary 
parts of (a(t)) Ta . We remark that result (351 agrees with an alternative derivation that does not use phase space 



techniques but employs the su(l,l) group structure of squeezing operators instead. 

Formula (351 demonstrates that the functional structure of a~ 1 {a(t)) TOC with respect to four initial state parameters 



{t, a} depends on just the three variables \a\, s, Atfi. At t = 0, the case Acf> — corresponds to phase squeezing: the 
uncertainty of the magnitude of (a) Ta is increased and that of its phase factor is reduced. An example for phase 
squeezing is the lighter ellipse shown in Fig. [3] On the other hand, A0 = ir corresponds to number squeezing: the 
uncertainty of |(a) TQ .|, which is the square root of the mean number of photons, is decreased. This is the case for the 
darker ellipse shown in Fig. [3] 

Note that solution & i(t\x) has different frequencies for different \x\. Thus the minimum uncertainty character of 
the initial state \ra) is lost for t ^ and restored at half-period multiples, £,LOqt = Nn. 

An important feature of the expectation value (351 is that for fixed s > it is a smooth bounded function in 
all variables; in particular it does not display the Moyal solution singularity when £u)2t approaches n/2. In order 
to interpret this we recall how the Heisenberg uncertainty principle works in the quantum phase formalism. Weyl 
symbols are often distributions and as such do not have any restrictions on localization or magnitude. For example, 
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a) 01 - b) 

FIG. 4: Mean a) position and b) momentum for number squeezed states (A0 = tt) for the case a = £ = 1. 




0.1 n 




a) 01 * b) 

FIG. 5: Mean a) position and b) momentum for phase squeezed states (A<j> = 0) for the case a = £ = 1. 




0.1 71 



the quantizer A(x') cf. (All is a bounded operator whose symbol is the delta function S(x' — x). The information 
about phase space uncertainty is encoded in the Wigner function [|ra)(ra|] (x). Only when the expectation value 
integral (29) is evaluated are the full effects of quantum uncertainty imposed. This phase space integration averages 
out the Moyal solution singularity giving a finite result. 

However, one may ask whether there is a surviving signature of the Moyal singularity at £w2^ = vr/2 in the 
observables (q(t)) Ta or (p(t)) Ta . For suitable values of r, a the answer is yes. In the next section we will present a 
general argument that these effects should be most evident for states similar to number-squeezed states, Acf> = tt. 
In Fig. |4] it can be seen that this is indeed the case: for strong number squeezing (s <C l,A</> = tt) position and 
momentum mean values are significantly enhanced around the singular point On the other hand, in the case of phase 
squeezing (A0 = 0) the singular behavior of the Moyal solution has virtually no effect on the mean amplitude (Fig.|5|. 
We remark that the peak in Fig. [5] appears at t = and therefore corresponds to the mean value of position and 
momentum in a squeezed coherent state in absence of the nonlinear interaction. 

Fig. [4] a) suggests that a singularity does appear in the limit of infinite number squeezing. However, for any fixed 
value £u>2t ^ tt/2 we have lim s _>o(ct(£))Ta = 0. On the other hand, for fixed s we have 



lim (a(t)) Ta = - exp 



This indicates that the peak becomes infinitely narrow in the limit of infinite squeezing and thus is not of physical 
significance. We remark that in current experiments a squeezing factor of about s = 0.1 can be achieved |25j . 

The general expectation value simplifies dramatically in several special cases. For purely coherent states, r = 0, 
expression (35 1 reduces to the well-known result (see, e.g., Eq. (40) of Ref. [2"rj] ) 



(&(t)>7 



a exp 



-iujit 



(36) 



This no-squeezing result does not show any evidence of the Moyal solution singulary. 
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Next consider the small £ behavior of (a(t)) Ta . In this limit all the noncommuting effects of the Heisenberg algebra 
for q,p are turned off. The Weyl symbol based semiclassical expansion is constructed by replacing 0oi with its 
semiclassical approximation ( 23 ) in the evaluation of the integral ( 34 ) . One finds 



(a(t)) Ta = a c i(t\x) [l + [2|T|e lA * - 2|a| : W(W + 4i|r| cos(A$))] £ + 0(£ 



(37) 



The leading £ — term term is just the complex statement of the Kerr classical flow (20) with the initial condition 



(q,p) = (q,p) = \/2 (Rea,Ima). The formula (37l may also be obtained by implementing a power scries expansion of 
(351. However the procedure using &Q{(t\x) has wider application in that it does not require an exact Moyal solution. 



The ^-linear correction factor senses the dependence on the squeezing |t| and the phase A<I>. For this expansion to 
be a good approximation to (a(t)) Ta the factor 1 in the curly bracket must be much larger than the correction terms. 
For \a\ = 1 this requires s S (0.8, 1.0) and \uj2t\ <C 1. Though this region of good approximation is an extremely small 
portion of the variable range shown in Figs. |4]and[5j it nevertheless covers a significant part of the experimentally 
accessible range. 



VI. FINITE EXPECTATION VALUES 



The result that the exact solution for the Weyl symbol of a well-established physical observable diverges periodically 
in time is surprising. Such time-periodic singularities are not present in the harmonic oscillator basis or in the 
classical solution for the Kerr Hamiltonian. In this section we therefore examine for general initial states whether 
such singularities could in principle survive the phase space averaging integral that defines an expectation value in the 
Weyl symbol picture. The results of the previous section show that the squeezed state expectation values of @oi(t\x) 
are finite for all times. Here we show that the expectation values of the general solution Q sm (t\x), whose amplitude 
diverges as (seci) s+m+1 , are finite for almost all possible quantum states. 

A very useful property of coherent states is their over-completeness, 



1 



1= — I \a)(a\d 2 a. 



We therefore can express the expectation value of an operator / as 

1 



Tr(p/) 



«) 2 



((3\p\a)(a\f\P)d z ad A (3 



with p the density matrix of the initial state of the system. Hence, to see if the singularities can appear for any 
quantum state it is sufficient to investigate the matrix element (a|/|/3) of an operator. Using Eq. (A3 I we can express 
this matrix element as 

(a\f\(3)= f f(x)(a\A(x)\f3)d 2 x, 



with / the Weyl symbol of /. In the following we will evaluate the integral over d 2 x for the set of Weyl symbols ( 16 I 
Using Eqs. (31 ) and (A2| it is not hard to see that, in complex coordinates z = q + ip, 



(38) 



where C = -[\a\ 2 + |/3| 2 + 2 i 5a*]/(2£). 

Using Eqs. (161,(38), d 2 x — \dzdz* and (A9) this leads to 



a\(a(tyy a(t) m \P) = I & sm (t\x) (a\ A(x) \(3) d 2 x 



— i(m— s)uj\ t 



(sec ty 



i+l i(2-s)t 



\ J (<x\&(x)\0) Jf Z ^j exp (-^*tan(i)^ dzdz* 



(sect) m+1 e l ^- s ^(a*) s i / (a\A(x)\p) exp f - l -zz* tan(£)J (-^=) dzdz*, 
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where we have performed a partial integration and used that (z* + £d z )(a\A(x)\{3) = y/2a* (a\A(x)\/3) . Performing 
the integration yields 

(a\(a{tyy a{t) m \/3) = (sect) m+1 e- i ( m - a ^ lt e c+i(2 - s)t " (^9 a ,\ J e -^ {1+ ' ltani)+ f {t3z ' +a ' z) dzdz* 

= (a*) S pm e -i(rn-s)t(u 1+ (rn+s-l)Su 2 ) exp J_ (j a |2 + ^2) + ^ e -3i(m-.)&*t^ _ 

This result is in perfect agreement with the corresponding expression derived from the solution of the ordinary 
Heisenberg equations of motion and thus demonstrates the consistency of the Moyal-Kerr solution (17 1. For our 
discussion it is important to note that none of these matrix elements contains a singularity. This means that the 
singularity can be considered as a feature of the Moyal representation of quantum mechanics that does not generate 
divergent expectation values. 

An intuitive explanation of why the singularity does not show up in expectation values is as follows. At the time 
t when the amplitude of Eq. (jlTj) diverges, the phase factor diverges as well. However, the diverging factor in the 
phase contains the mean photon number zz* . Consequently, the phase divergence is different for states with different 
photon numbers. Therefore, for t sufficiently close to the singular value tt/2, the phase factor would average to zero 
for any state that has a variance in the number of photons. 

Hence, the only chance to see the singularity would be in a state where the number of photons is exactly, known, i.e., 
number states. But number states correspond to states for which the phase is completely undetermined, so that any 
expectation value with s =/= m would average to zero for any value of t. Hence, the uncertainty relation An A<p > 1/2 
(see, e.g., Ref. [27 ) for photon number and phase prohibits the appearance of the Moyal singularity in a quantum 



photon observation (17 1 



VII. CONCLUSION 



In this paper we have derived an exact solution (16 1 for the Weyl symbol phase space representation of the Kerr 
model of nonlinear quantum optics. This solution exhibits singularities that are absent for the classical solution and 
depend on the intensity of the light field. This intensity-dependence guarantees that expectation values for all initial 
states remain finite. On the other hand, a signature of the singularity appears in the form of finite peaks in the 



expectation values (35) for number-squeezed states of light. A number of open questions remain. The single- mode 
Kerr model that we have studied gives only a good description for photons in optical cavities of extremely high 
finesse. A generalization of our results for a multi-mode theory of propagating photons would therefore be desirable. 
Alternatively, an imperfect cavity could be modeled by studying a Kerr model that is coupled to the environment and 
exhibits Langevin noise. Both of these aspects have been addressed by Kartner et al. [25] in the context of a specific 
noise model using the Wigner function. The corresponding Weyl symbol representation of the Heisenberg-Langevin 
equations of motion may be free of divergences and would open new perspectives for phase space descriptions of 
quantum systems. 
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APPENDIX A: WEYL SYMBOL QUANTUM MECHANICS 



This appendix summarizes properties of the quantum phase space method that are employed in this paper. We 
collect the various Weyl symbol identities in a notation suitable for quantum optics. The account below closely 
matches that found in Refs. fl9| [29] . 

For a single-mode photon and a suitable choice of reference point, the electric field strength E is proportional to 
a + a) (or equivalently to q) and the magnetic field strength B to i(a — SJ) or p. For this reason the phase space base 
manifold for a single mode state is the real line R. The noncommutivity of q and p arise from the mode operators a 
and aJ . The quantum state space is that spanned by the harmonic oscillator basis, or equivalently the Hilbert space 
of one dimensional square integrable wave functions, 7i — L 2 (M.,C). Likewise, the associated classical phase space is 
T*M. — M 2 equipped with the standard Poisson bracket. 
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Weyl quantization maps functions on T*M. into operators on Ti. A unified characterization of both quantization 
and de-quantization is achieved via a quantizer [30 1 3T] 132] . Let {A(x) : x = (q,p) € T*K} be a x-depcndcnt family 
of bounded, self-adjoint operators on Ti defined by their action on a wave function, ip 

^V) = s ^ ex P - g)) v(2g - g') , (Ai) 

or, equivalently, as an integral kernel 

(g'|A»|g") = i exp Qptf - q")\ S(2q - q' - q") . (A2) 

Then both quantization and dequantization are constructed from A(g,p) via 

/= I f(x)A(x)d 2 x, [/] w ( I ) = (2<)Tr/A(x) (A3) 



where Tr is the trace on Ti. This pair of linear transformations are mut ual inverses, so / = [/] w . The notation [/] w 
indicates the Weyl symbol of the operator /. The second identity in Eq. (A3 1 is proportional to the Wigner transform 
of/, c/. §. 

This bijective correspondence between phase space functions and operators is simple in a variety of important cases. 
For example, operators /(<?), g(p) and the identity on Ti have symbols f(q), g(p) and the constant function 1. The 
quantizer has symbol, [A(x)] w (x') — S(x — x'); in turn this implies that the exponential operator e lu ' x , u £ M 2 has 
the symbol e m ' x . 

The Weyl symbol framework is a Hilbert algebra, namely a complete linear space L with three basic structures: an 
associative product *, an involution * and an inner product The product of operators on Ti is mirrored by the 

noncommutative product of Weyl symbols. This star product is defined by / * g = [f g] w - Given the * product, the 
Moyal bracket is defined as 

{/, 9}m = ~ [/, g]w = ^ (/ ★ g - g * f) ■ (A4) 

The ★ product has three useful representations. 
The first is Berezin's integral form [3"3"] 

f*g(x) = ^j^ 2 J J f(x 1 )g(x 2 )exp^-(x 1 Ax 2 + x 2 /\x + xAx 1 )^d 2 x 1 d 2 X2. (A5) 

Here x\ A x 2 = X\ • Jx 2 . 

Next is Groenewold's derivative expansion [34] 

/ * g(x) = exp (^di ■ Jd 2 ^j <f,g> (x) . (A6) 

Above di and d 2 are gradients acting on the first (/) and second (g) arguments of the product -<f,g >■ The Poisson 
bracket {/, g}, in this notation, is d\-Jd 2 <f,g >~ followed by diagonal evaluation, x\ = x 2 — x. In the case where 
/, g are £ independent and suitably smooth, the series expansion of Eq. (A6l gives a ^-asymptotic expansion of f * g 
with the Poisson bracket term as the leading semiclassical correction. 
The third realization of the * product is the left, right form. Define 

L = x+^Jd x K = x-^Jd x , (A7) 

then for smooth /, g 

f * g (x) = (/(L) g) (x) = {g(K) f) (x) . (A8) 

The differential operators C and 1Z commute. 

The involution operation on L is complex conjugation. It is the symbol analog of the adjoint operation on Ti. If 
/ = [/] w then /* = [/t] w . 
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For suitably restricted f,g (e.g. both Hilbcrt Schmidt), the trace of the operator product defines the L inner 
product, in detail 

(2^)Tr/5= / f*g(x)d 2 x = f f(x)g(x)d 2 x = (f\g) L . (A9) 
Jr 2 Jr 2 

This formula shows that quantum expectation values in L are obtained by phase space integration. For example, this 
occurs if / is a density matrix and g is any observable. 

Weyl quantization has a simple covariance property. A unitary operator V is called metaplectic if VxV' — Sx for 
some symplectic matrix S, i.e. SJS T = J. If / has the symbol /, the affine canonical covariance property [TS] is the 
statement that 

[VfV^ w (x) = f(Sx). (A10) 
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